Vortex rectification effects in films with periodic asymmetric pinning 
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We study the transport of vortices excited by an ac current in an Al film with an array of 
nanoengineered asymmetric antidots. The vortex response to the ac current is investigated by 
detailed measurements of the voltage output as a function of ac current amplitude, magnetic field 
and temperature. The measurements revealed pronounced voltage rectification effects which are 
mainly characterized by the two critical depinning forces of the asymmetric potential. The shape 
of the net dc voltage as a function of the excitation amplitude indicates that our vortex ratchet 
behaves in a way very different from standard overdamped models. Rather, as demonstrated by the 
observed output signal, the repinning force, necessary to stop vortex motion, is considerably smaller 
than the depinning force, resembling the behavior of the so-called inertia ratchets. Calculations 
based on an underdamped ratchet model provide a very good fit to the experimental data. 

PACS numbers: 05.40.-a, 74.78. Na., 74.40.+k, 85.25.-j 



From the point of view of classical thermodynamics, it 
is not possible to induce directed motion of particles by 
using equilibrium fluctuations only, otherwise it would 
constitute a perpetuum mobile of the second kind 0. 
Nevertheless, non-equilibrium fluctuations, such as pe- 
riodic excitations or a "colored" noise, are allowed to 
take advantage of the asymmetry of a periodic ratchet 
potential to promote motion of particles in a preferen- 
tial direction 0|. New solid-state-based ratchet systems 
are currently being developed for controlling the motion 
of electrons and fluxons, as well as for particle sepa- 
ration Q and electrophoresis In particular, ratchet 
potentials in superconducting devices may be very useful 
to control the dissipative motion of fluxons, which causes 
undesired internal noise. 

Modern lithographic technics make it possible to fab- 
ricate periodic arrays of vortex pinning sites with size 
and shape that can be easily tuned, thus giving an inter- 
esting perspective for making different asymmetric pin- 
ning potentials. In this context, several ideas to control 
flux motion by applying an ac excitation have been pro- 
posed [a S H @ > but up to now only a few experiments 
have been realized [lCJ, . One realization has been re- 
cently implemented on a Nb film with a square array of 
nanoscopic triangular magnetic dots 0. The authors 
reported rectification of the ac driven vortices due to the 
asymmetric shape of the dots. Nevertheless, the detailed 
dynamics of vortices in such structures is not yet com- 
pletely understood. 

In this Letter we investigate a composite square array 
of pinning sites, with its unit cell consisting of a small 
and a big square antidot separated by a narrow super- 
conducting wall, as a vortex rectifier. As demonstrated 
by our dc and ac transport measurements at several fields 
and temperatures, this configuration is able to break the 



reflection symmetry of the total effective pinning poten- 
tial and promote flux quanta rectification. Moreover, 
our data reveals a remarkable hysteresis in the current- 
induced pinning-depinning process. This gives an appar- 
ent inertia to the system with important consequences in 
the overall dynamics. Finally, we discuss how standard 
overdamped models fail in describing these results and 
propose an underdamped ratchet model that provides a 
very good fit to our experimental data. 

The pattern was prepared by electron-beam lithogra- 
phy on a SiC>2 substrate. The superconducting thin film 
is a 38 nm thick Al, in a 5 x 5 mm 2 cross-shaped geome- 
try [see Fig. 0(a)] to allow four-point electrical transport 
measurements in two perpendicular current directions. 
The central part of the cross consists of two 300-^m- 
wide strips containing the nanoengineered array (period 
a p = 1.5 |xm) of asymmetric pinning sites [see Fig. 0(b)]- 
This gives a value of 0.92 mT for the first matching field, 
Hi = §o/ a p (here $0 = h/2e is the flux quantum). In 
each strip, the voltage contacts are 2 mm apart. The 




FIG. 1: Layout of the Al film, (a) Cross-shaped geometry of 
the sample to allow for transport measurements in the x and 
y directions, (b) atomic force micrograph of a 5 x 5 )im 2 area 
of the asymmetric pinning sites 
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FIG. 2: (a) Contour plot of the magnetic field and dc current 
dependence of the voltage difference for dc currents applied 
in the x and -x direction V(I- X ,H) - V(I X ,H) at T/T c = 
0.973 (b) Contourplot of the net dc voltage Vda(Iac,H/Hi) 
as a function of the magnetic field and amplitude of the ac 
sinusoidal current at a frequency 1 kHz and T/T c = 0.973. 



sizes of the small and big antidots are 300 x 300 nm 2 and 
600 x 600 ran 2 with a separation between them of 90 nm. 
The magnetic field was applied perpendicularly to the 
sample surface. The superconducting critical tempera- 
ture (T c = 1.469 K, AT = 3mK) was obtained by using a 
resistance criterion of 10% of the normal state resistance. 
From the residual resistivity at 4.2 K, the elastic mean 
free path (l e i = 3.9 nm) was found. From this we can 
calculate the coherence length £(0) = 67.5 nm and the 
penetration length A(0) = 195 nm using the dirty limit 
expressions. We can therefore conclude that the pat- 
terned film with the composite array of square antidots is 
a type-II superconductor (k = = 2.89 > -j= = 0.71). 

Fig-EI a ) shows a contour plot of the magnetic field and 
current dependence of the difference in induced voltage 
for the dc currents flowing in the positive and negative 
x direction, V(I X ,H) - V(I- X , H), at 97.3 % of the crit- 
ical temperature. The voltage drop in the x direction 
is a measure of the vortex velocity in the y direction. 
Since the I_ x and I x dc currents induce a Lorentz force 
in the negative and positive y directions, the vortices will 
probe the asymmetry of the pinning sites. We see that 
for all fields (0 — l.lf/i) there is a dc voltage difference 
confirming the presence of an asymmetry in the pinning 
potential. This difference is the largest for fields just be- 
low the first matching field, Hi, which can be explained 



by the high symmetry of the vortex configuration which 
cancels out the vortex- vortex interactions. In this case 
the resulting forces acting on vortices are coming only 
from the asymmetric pinning potential. To get the data 
shown in Fig.^b), a 1 kHz sinusoidal current was applied 
to the sample and the output signal was measured with 
an oscilloscope used as an analog-digital converter. Inte- 
grating this signal gives the average dc response, which 
is plotted as a function of the normalized magnetic field, 
H/Hi, and amplitude of the applied current, I ac . Again, 
we can clearly see a ratchet effect at all fields from 
to 1.1H\. This effect is maximum just below H\ and 
i-ffi, where vortex- vortex interactions are also mostly 
cancelled. The measurements carried out at other tem- 
peratures (T/T c = 0.967 up to T/T c = 0.98) show similar 
features. Nevertheless, the rectification effect is weaker 
at higher temperatures, which can be explained by an 
increase in £(T) and A(T), making the size of the vor- 
tices big compared to the antidot size. In this case the 
vortices cannot sense the asymmetry anymore. We have 
also repeated the measurements for several frequencies in 
the range 10 Hz - 10 kHz. No significant change in the 
results was observed, which means that our experiments 
are carried out deep in the adiabatic regime [2|. 

As it is clearly seen from the vertical scans in Fig. 0b) , 
the net dc voltage increases sharply with the ac-current 
amplitude up to a maximum value and then decays 
smoothly to zero. This is a typical behavior of ac-driven 
objects in a ratchet potential in the adiabatic regime. 
At low frequencies, the net dc flow (v) of particles in- 
creases monotonically when the excitation amplitude A 
is between the weaker, fdi, and the stronger, fd2, depin- 
ning forces of the asymmetric potential. We shall refer 
to this region as the rectification window. Since in this 
amplitude range the force cannot overcome the potential 
barrier for negative force direction, motion occurs only 
during the positive half loops of the ac excitation, that 
is, the system behaves as a half-wave rectifier. For am- 
plitudes A > fd2, the particles are allowed to travel back 
and forth, thus leading to a decrease in the rectification 
signal, which vanishes slowly at high drives as the ratchet 
potential is gradually averaged out by the fast-moving 
particles. It is worth noting that these general properties 
of the ratchet effect also depend on the shape of the ac 
excitation 0- For instance, square- wave excitations lead 
to sharper rectification effect and a shorter tail in the (v)- 
A characteristics as compared to sinusoidal excitations. 
Since the regime of the dc measurements is equivalent to 
irradiating the system with a square wave pulse, the dc 
data presented in Fig. HJa) looks much sharper than the 
ac data shown in Fig.[5fb). 

A closer look at the ac-amplitude dependence of the 
net dc voltage reveals additional and important features 
of the ratchet dynamics in our system. Figs.|21and0]show 
detailed dc voltage versus ac amplitude (Vd c -I a c) charac- 
teristics. For fields below the first matching field some 
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FIG. 3: Net dc voltage as a function of ac amplitude at 
H/Hi = 0.96 and T/T c = 0.973. The current is normal- 
ized by the Idi = 473 /J.A. Inset: Voltage output at the ac 
amplitudes (I) 477/iA, (H) 502 fiA, (III) 544 and (IV) 
704 piA; (a) time evolution of the vortex velocity v(t) from an 
overdamped molecular dynamics simulation (see text); and 
(b) v(t) calculated from the inertia ratchet model. The hor- 
izontal dashed lines indicate the instantaneous value of ex- 
citation where vortices are rfepined [F(t) = Fd] and repined 
[F(t) = F r ]. The subscripts 1 and 2 stand for the weaker and 
the stronger depinning forces respectively. 



Vdc-Iac curves, as the one at H = 0.96i/ c i (Fig. |3J), have 
two peaks. Detailed measurements of the time evolution 
of the output voltage V(t), indicated by roman numbers, 
suggest a superposition of two ratchet effects. Rectifica- 
tion is first triggered just above I ac = 473 /zA and the 
output signal is present only at the positive half-loop of 
the ac current, as shown in panel I. At I ac — 495 /zA, a 
stronger ratchet is triggered. It can be identified in panel 
II as a second jump in the signal in the positive half-loop. 
Right above I ac = 520 /zA a small signal in the negative 
half loop appears (sec panel III), thus demonstrating that 
the weaker ratchet is already outside its rectification win- 
dow. In panel IV, the signal is strong in both directions, 
indicating that the driven vortex lattice moves back and 
forth as a whole. A possible reason for the presence of two 
ratchets can be a complicated plastic dynamics in which 
vortex channels flow along specific rows of pinning cen- 
ters for currents considerably lower than the depinning 
current of the whole vortex lattice. Plastic dynamics was 
thoroughly studied by numerical simulations in periodic 
arrays of symmetrical pinning centers The local 

decrease in pinning efficiency in a given pinning row is 
caused by discommensurations (vacant pinning sites for 
fields H < Hi) distributed along the row. The extension 
of these ideas to arrays of asymmetric pinning sites seems 
to be rather straightforward. At fields very close or at Hi 



this effect is minimized. At H = 0.98-ffi for instance, the 
Vdc-Iac characteristic, shown in Fig. 0] presents a single 
well-defined rectification peak. 

Another remarkable feature disclosed in the V(i) data 
is that the repinning force f r necessary to stop vortex 
motion in either direction is always smaller than the cor- 
responding depinning force fd [see Fig-El insets (I)-(IV)]. 
In other words, the V-I characteristics is hysteretic. This 
is a very robust phenomenon which has been observed in 
the whole field, frequency and temperature ranges used 
in our experiment and, to the best of our knowledge, has 
not been reported before. To understand the experimen- 
tal results, we have performed molecular dynamics (MD) 
simulations of the overdamped equations of motion for 
vortices in a square array of asymmetric pinning poten- 
tials at fields close to the Hi and subjected to a sinusoidal 
Lorentz force of very low frequency (the details of the 
model and simulations are reported in Ref. [t|). In the 
whole range of excitation amplitudes, no hysteresis was 
observed (/,■ = fd), the v(t) curves at each half period 
are symmetrical, in contrast with the clearly asymmet- 
rical shape of the V(t) experimental data. In addition, 
the (v)-A characteristics look qualitatively different from 
the experimental results. A possible reason for the mis- 
fit between the MD calculations and the experimental 
results is that the depinning process of an overdamped 
particle is reversible because it carries no significative in- 
ertia and cannot be deformed. However, vortices are soft 
objects that can be deformed under the action of com- 
peting forces and they also carry a small mass. In the 
dirty limit, the vortex mass can be several times smaller 
than the electron mass. In this case, inertia is only rel- 
evant at frequencies typically above 10 12 Hz Since 
our sample is in the dirty limit, we believe the contri- 
bution of vortex mass is negligible in our experiment. 
Therefore, the apparent inertia is more likely due to vor- 
tex deformation. Indeed, recent numerical simulations of 
the Ginzburg-Landau equations 0] suggested that the 
current-induced depinning of vortices in a periodic array 
of antidots is dominated by strong elongation of the vor- 
tex cores towards a neighboring antidot. Interestingly, 
the authors observed hysteresis in the pinning-depinning 
process resulting from this behavior. 

To illustrate how this hysteresis influences the overall 
ratchet dynamics, we model vortices still as rigid parti- 
cles but carrying an apparent mass M. For simplicity, 
we consider that the vortex lattice is perfectly commen- 
surate with the pinning array and therefore we shall ana- 
lyze the experimental results measured very close to the 
first matching field (H = 0.978i?i) shown in Fig. H In 
this case, due to the cancellation of the vortex-vortex in- 
teractions, the vortex dynamics can be reduced to that of 
one particle moving in a one-dimensional potential. The 
corresponding equation of motion is 
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FIG. 4: Normalized dc voltage as a function of ac am- 
plitude. Circles: experimental data at H/Hi = 0.98 and 
T/T c — 0.973. Line: numerical fit using the inertia ratchet 
model [Eq. Q] with only one fitting parameter, 7 = f .8. In- 
set: numerical integration of the inertia ratchet model for dif- 
ferent 7 (symbols) and the overdamped limit solution (line). 

where rj is the viscous drag coefficient, v is the instanta- 
neous vortex velocity and U p (x) is the one-dimensional 
pinning potential, which we choose in a simple form as 
U p {x) — — {7[sin(27rx/a p )+0.025sin(47rx/a p )]. This gives 
an asymmetry ratio fdi/ fd2 = 0.925 similar to that ob- 
served in our data shown in Fig. 0] (Idi — 480 fj,A and 
Id2 = 519 /LtA). In order to be closer to the experimen- 
tal conditions, we will restrict ourselves to the adiabatic 
regime, that is, we shall consider only frequencies much 
lower than the libration frequency of the pinning poten- 
tial uj p ~ 2ir^Ju /Map and the viscous drag to mass ratio 

7 = -q/M. In this way, by normalizing the ac amplitude 
by the first depinning force and the dc response by max- 
imum response, the only fitting parameter for comparing 
the calculations with the experimental data is 7. 

The overdamped limit corresponds to M — > 0, or, 
equivalently, to 2-7T7 3> w p . In this limit a semi-analytical 
approach may be used. The average flow of vortices 
in one cycle of the ac excitation is given by (v) = 
■p f Q dtv(t) where P = lo/2tt is the forcing period and 

«(*) = nf a, ( 2 ) 
'JO -U^(x)+Asin(ujt) 

This model reproduces exactly the MD simulations for 
H = H\. The j{v)-A characteristic calculated for the 
overdamped limit is given by the full line in the inset of 
Fig. 0] The curve presents a long and slowly decaying 
tail qualitatively different from the curves observed in 
our experiments. 

For finite 7 values, we integrate Eq.n numcr i ca lly by a 
finite difference method. The results for several 7 values 
(in units of uj p /2it) are shown in the inset of Fig.0] For 



7 = 4, the ry{v)-A characteristics is already very close to 
the overdamped limit. For smaller 7 values, the effect 
of inertia becomes increasingly more pronounced. The 
curves look sharper and more restricted to the rectifica- 
tion window. As shown in the main frame of this figure, 
the 7 = 1.8 curve provides an excellent fitting to the ex- 
perimental data. The time evolution of the velocity v(t) 
for 7 = 1.8 is shown in the inset (b) of Fig. Olfor a few 
ac drive values. Note that the curves correctly reproduce 
the asymmetry observed in the experimental V(t) data 
[insets (I)-(IV)], with F r ^ F d . The v(t) curves calcu- 
lated for 7 > 4 resemble the overdamped limit [Eq. , 
which gives exactly Fd — F r and no hysteresis. 

In conclusion, we have demonstrated voltage rectifi- 
cation in a superconducting film with a composite peri- 
odic array of asymmetric pinning sites. The rectification 
was observed in dc and ac measurements in a broad field 
range and for different temperatures, showing that our 
sample provides a two-dimensional asymmetric potential 
which is able to control vortex motion. For fields below 
the first matching field, our data revealed a competition 
between two ratchet dynamics, which can be understood 
as the competition between plastic channelling dynam- 
ics and coherent vortex motion. We also have observed 
a strong irreversibility in the pinning-depinning process, 
which we believe is due to current-induced elongation of 
the vortex core. Calculations based on an underdamped 
ratchet model provided a very good fit to the experimen- 
tal data and showed that the hysteresis found in our ex- 
periment gives rise to an apparent inertia with important 
consequences to the ratchet dynamics. 
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